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We have theoretically investigated the electromechanical properties of a freely suspended carbon 
nanotube that is connected to a constant-current source and subjected to an external magnetic field. 
We show that self-excitation of mechanical vibrations of the nanotube can occur if the magnetic 
field H exceeds a dissipation-dependent critical value which we find to be of the order of 10-100 
mT for realistic parameters. The instability develops into a stationary regime characterized by time 
periodic oscillations in the fundamental bending mode amplitude. We find that for nanotubes with 
large quality factors and a magnetic- field strength just above He the frequency of the stationary 
vibrations is very close to the eigenfrequency of the fundamental mode. We also demonstrate that the 
magnetic field dependence of the time averaged voltage drop across the nanotube has a singularity 
ai H = He. We discuss the possibility of using this phenomenon for the detection of nanotube 
vibrations. 



Nanoelectromechanical systems (NEMS) containing a 
suspended carbon nanotube (CNT) vibrating at radio 
frequencies (RF) have received increasing attention re- 
cently. Advantages of CNT mechanical resonators in- 
clude their high resonance frequencies — up to the GHz 
range — - their low dissipative losses and the possibility 
to tune the resonance frequency by adjusting the tension 
in the tube^'^. CNT based NEMS devices have already 
shown a great potential for a plethora of technological 
applications including mass sensing ' and tunable high 
frequency electronics ' ' ' . However, most of the devices 
which have been realized thus far are passive resonators 
which perform frequency filtering of the incoming RF- 
signaP''^'^'^. Here we propose an active oscillator based 
on an current biased doubly-clamped CNT for which the 
conductance depends monotonically on the nanotube de- 
flection. Such a deflection sensitive resistance has been 
demonstrated for a semiconducting single-walled CNT 
suspended over a gate electrode^''^'^. The active feedback 
is provided by a Lorentz force induced by a constant mag- 
netic field directed perpendicular to the direction of the 
current and parallel to the gate electrode. We show that 
, by applying a constant external current in a sufficiently 
high magnetic field, we obtain mechanical instability and 
self-sustained mechanical oscillations at a frequency close 
to the mechanical resonance frequency. Furthermore, we 
show that the mechanical instability results in oscillations 
of the voltage drop across the nanotube accompanied by 
a deviation of the time averaged voltage from its static 
time-independent value. 

The proposed oscillator device is shown in Figure 1. 
A semiconducting carbon nanotube is suspended over a 
gate electrode and subject to a homogeneous magnetic 
field H perpendicular to the nanotube and parallel to the 
gate electrode. The gate-nanotube separation depends 
on the deflection z{t^x) of the tube from its equilibrium 
straight configuration {x is the coordinate along the nan- 
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FIG. 1: Sketch of the proposed active oscillator device. A 
semiconducting carbon nanotube is suspended over a gate 
electrode and connected to an external dc current source. An 
external magnetic field, applied perpendicular to the direction 
of the current, gives rise to a Lorentz force that deflects the 
tube towards the gate, which affects the resistance and pro- 
vides a feedback mechanism that for large enough magnetic 
fields leads to self-sustained nanotube oscillations (see text). 
The inset shows an equivalent electric circuit 



otube). We will assume that the nanotube mechanics 
is completely characterized by the amplitude of the fun- 
damental bending mode u{t)^ and consequently let the 
time dependence of the mechanical deflection have the 
form z{t^x) = u{t)ipo{x)^ where (po{x) is the normalized 
profile of the fundamental mode . This assumption cap- 
tures the essence of the physics behind the phenomena 
to be considered below, for which the dependence of the 
mutual capacitance Cg{u) between nanotube and gate 
on the amplitude of the fundamental mode is central. 
This is because the concentration of charge carriers in 
the semiconducting nanotube p{CG{u)^Vg) is a function 
of both the gate voltage Vg and the mutual gate capaci- 
tance Cg{u)^ which means that ultimately the nanotube 
resistance R{p{Cg{u)^ ^g)) will depend on the nanotube 
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deflection. In the experiments of Sazonova et al^, this ef- 
fect was utihzed for the detection of resonant mechanical 
excitations of the nanotube. 

When an external current-source feeds a current Iq to 
one of the leads, a Lorentz force proportional to the cur- 
rent through the nanotube and the magnetic field causes 
a deflection of the CNT. For displacements on a scale 
much smaller than the length of the wire, its curvature is 
negligible and consequently one can consider the Lorentz 
force to be the same at every point. Since the force is 
almost uniform over the wire, it is mainly the fundamen- 
tal bending mode that will be affected, which provides 
another justiflcation for our assumption that all other 
modes can be neglected. Assuming the nanowire to be 
an elastic beam whose motion can be described by lin- 
ear continuum mechanics, we arrive at the following set 
of equations for the time evolution of the voltage drop 
V{t) across the nanotube and the amplitude of the fun- 
damental mode u{t): 



mu(t) + -fiiit) + nuit) = LHI 
CV{t) =Io-I, 1 = V/R{u(t)) , 



(1) 



Here m and L are the effective mass and length of the 
suspended part of the nanotube, n and 7 are the effective 
spring and damping constants, and C is the capacitance 
of the junction. 

The system of equations (1) has a time independent 
solution given by 



u{t) = uq = LHIq/ n 

V{t) = Vo = R [^'^ lo. 



(2) 
(3) 



A linear stability analysis of this solution yields a secular 
equation for the Lyapunov exponents A of the form 

P(A,/3) = (A2 + g-ic^oA + c^g)(A + c^^)-/3c^g^ = 0, (4) 



where co'o = ^J Kjm is the eigenfrequency of the funda- 
mental mode, ujR = 1/R{uq)C^ and Q = is the 
quality factor characterizing the damping of the mechan- 
ical vibrations due to its interaction with the thermody- 
namic environment. The parameter 
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is proportional to the Lorentz force and will be referred to 
as the electromechanical coupling parameter of the sys- 
tem. It is worth to note that since uq^ and consequently 
Ir^ is a function of the magnetic fleld and the external 
current, /3 is not necessarily directly proportional to the 
externally adjustable parameters H or Iq. How much 
this affects the possibility to control the electromechan- 
ical coupling through the magnetic fleld depends on the 
speciflc geometry and the values of other parameters such 
8isVg. It turns out that there is a simple relation between 
the value of the coupling parameter /3 and the sign of the 



differential resistance in the static regime. Indeed, from 
the time independent solution (3) we obtain the relation 
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from which we see that the differential resistance in the 
time independent regime is positive if /3 < 1, while for 
/3 > 1 it is negative. 

For a high quality factor, Q ^ 1, and weak electrome- 
chanical coupling, /3 <C 1, the (approximate) solutions of 
equation (4) — together with a third root which is always 
real and negative — are 



" (7) 

From equation (7) follows that the two Lyapunov expo- 
nents Ai^2 have a positive real part when 
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(8) 



This means that for large quality factors, Q ^ 
max(co'o/cc;i?, cji^/cjo), the regime of time independent 
charge transport through a static nanotube becomes un- 
stable with respect to self-excitation of mechanical vi- 
brations, if the coupling parameter j3 exceeds the critical 
value j3c' The natural way to reach this self-excitation 
regime is to increase the magnetic fleld since increas- 
ing /o could lead to overheating. Hence we choose to 
characterize a particular setup by the critical magnetic 
fleld He above which the system will be unstable. From 
equation (8) follows that for large quality factors and 
'^0 = ^R the critical magnetic fleld is 

He = 2i^iR/QLIo. 

If we estimate the characteristic length iji as the distance 
between nanotube and gate electrode and take £r ^ 0.1- 
1 /j.m, 1 /j.m, K.^ 10-^ N/m, Q ^ 1000 and Iq ^ 0.1 
/iA, we obtain a critical magnetic fleld of the order He ~ 
10-100 mT. 

In the above section we analyzed the situation when 
the mechanical friction is very small. It turns out that 
the time independent regime becomes unstable for large 
enough electromechanical coupling even if the mechani- 
cal subsystem is overdamped, that is if Q <C 1. Indeed, 
the system becomes unstable when the real part of the 
Lyapunov exponents changes sign from negative to pos- 
itive. It means that at a certain critical point j3 = j3c 
the secular equation (4) has two purely imaginary roots 
Ai^2 = ^ioOc These are obtained from the two equations 

ReP(iUc, Pc) = ^R^o - ^li^R + Q"^^o) = (9) 
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whose general solution is 

1 UJqUJr + (5(cJo +^|) 



QujR 



(10) 



Important to note is that pc is always positive. In partic- 
ular, from this follows that the instability can only occur 
when the direction of the magnetic field is such that it 
deflects the nanotube towards decreasing resistance. In 
our geometry it means that the Lorentz force must be 
directed towards the gate electrode. 

From equations (10) it follows that for small quality 
factors, Q <C max(l, cjo/cj^^), and when the coupling 
parameter j3 has just overcome the critical value j3c = 
^r/{Q^o)i the static regime becomes unstable with re- 
spect to oscillations of angular frequency ujc ^ VojqujrQ^ 
which is always smaller than the mechanical resonance 
frequency. It is worth noting that, for an overdamped 
mechanical system the critical value of the coupling pa- 
rameter is bigger than unity, and as discussed previously 
this implies that the instability in this regime could be 
considered to be caused by a negative differential resis- 
tance. Estimating the critical magnetic field with the 
same parameters as before, but now with Q = 0.1, we 
obtain He ~ 100-1000 T, which is too large for exper- 
imental observation. Therefore, for the remaining part 
we will focus our attention on the low-dissipation regime, 
which can be achieved at least at low temperatures ' . 

The way in which the instability evolves is largely de- 
pendent on the magnitude of /3. For a coupling pa- 
rameter just above the threshold, Pc ^ P <C 1, we 
may analyze the development of the instability with the 
Ansatz u{t) = uq -\- A{t) sm{ujot)^ assuming A{t) to be 
a slowly varying function on the scale 1/uJo- Substitut- 
ing this Anzats into Eq. (1) and averaging over the fast 
oscillations , we obtain an equation for A, 



A{t) = aiuJoA{t) 



(11) 
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From equation (11) follows that for j3 > j3c {H > He) 
an initially small amplitude A(t = 0) <C 2£r will at first 
increase exponentially in time. There are then two differ- 
ent scenarios for the further development, depending on 
the sign of the coefficient bi. If 6i < and /3 — Pc 




0.02 0.04 0.06 0.08 0.10 

(J3 -MPc 



FIG. 2: The above graph show the results of a number of 
computer simulations of the system of equations (1) aimed 
at illustrating the regimes of soft and hard instability. In all 
the simulations we chose the quality factor Q = 10^ and the 
resistance R{u)/Ro = (1 + e"^^^"^o^/^^)/2. The stationary 
mechanical deflection is then equal to the coupling parame- 
ter: uo/in = /3. Initial values were u{0) = u{0) = V{0) = 0. 
In the first series of simulations we put uor = cjq for which 
formula (12) predict soft instability, the solid circles mark the 
saturation amplitude of the mechanical oscillation for different 
values of the ratio {/3 — pc)/ Pc- The results in this first series 
of simulations are in very good agreement with the analytical 
estimates. In the second series of simulations, marked with 
the solid squares in the graph, we chose ujr = 2(jJo for which 
formula (12) predicts hard instability. In this regime the sat- 
uration amplitude appears to have a discontinuous jump at 
the threshold of instability. Two curves are fitted to the data 
points. 



the amplitude saturates at the value 



As = 21 R. 
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(13) 



which vanishes for /3 = /3c {H = He) and corresponds to a 
"soft" instability. For hi > the amplitude may saturate 
at a value, which does not vanish as /3 ^ /3c (from above). 
This is called a "hard" instability. The different scenarios 
of soft and hard instabilities are illustrated in Figure (2). 

We note in passing that in the case of a hard instabil- 
ity, self-excitation results in a large amplitude ~ Ir of 
the mechanical oscillations and hence nonlinear effects in 
the nanotube dynamics become important. Therefore, in 
order to analyze this regime in detail one would have to 
consider nonlinear elastic mechanics. This would, how- 
ever, take us too far outside the scope of the present 
work. 

The result of the instability discussed above is a sta- 
tionary regime that is not only characterized by an os- 
cillating mechanical displacement of the nanotube but 
also by an oscillating voltage across the trench. In the 
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FIG. 3: The above graphs show the result of a computer 
simulation of the system of equations (1) aimed at illustrating 
the phenomenon of the voltage drop. For this purpose we 
chose the quality factor Q = 10^, the RC-frequency ur = cjq, 
the coupling parameter /3 = 1.1 /3c = 0.0219 and the resistance 
R{u)/Ro = (1 + e-2(--^o)/^i?)/2, where m/in = f3. Initial 
values were u{0) = u{0) = V{0) = 0. As can be clearly seen, 
the time averaged voltage deviates more and more from Vb 
as the amplitude of the mechanical oscillation increases. The 
actual oscillations are so fast on the timescale used that they 
are not resolved in the graph. 

regime of high quality factors, or more precisely when 
ujR ^ C(Jo/Q, the oscillation frequency will be very close 
to the mechanical resonance frequency a;o/27r, typically 
of the order 0.1-1 GHz. Detecting such fast voltage oscil- 
lations may pose serious experimental challenges. The 
nonlinear ity of the system does, however, result in a 
nonzero deviation of the time-averaged voltage V{H) 
from the static value Vo{H). This phenomenon is illus- 
trated graphically in Figure (3). If Ag <C in the time- 
averaged voltage can be approximated to lowest order 
with the formula 

V{H) - VojH) _ fdin 

Vo{H) \duo^ cv%+cv^J \2eRj ^ ' 

Thus, in the case of a soft instability the saturation 
amplitude can be calculated from formula (13) and is 
found to be proportional to {H — He) /He. Accordingly, 
in this case we have the following relation: 

^^o^'^eiH-H^) (15) 
where 0{x) is the Heaviside step function. Thus, at the 



critical point H = He 8i soft instability manifests itself as 
a jump in the derivative dV/dH whereas a hard instabil- 
ity would show up as a discontinuous jump in V. Hence, 
the instability could be detected by measuring the av- 
erage voltage drop across the suspended nanotube as a 
function of magnetic field. 

Important to note is that for our geometry the charac- 
teristic length is of the order of the distance between the 
nanotube and the gate. Therefore one can expect strong 
nonlinear effects in the nanotube dynamics to start to 
dominate at an amplitude which is just a small fraction 
of the characteristic length. Hence the validity of formula 
(11) is rather limited. However, we believe that formula 
(14) could still be valid as an approximate relation be- 
tween the voltage drop and the saturation amplitude. At 
the other extreme one could consider systems with a very 
short characteristic length, for example with an STM-tip 
positioned above the CNT (see for example Le Roy et 

14,24^^ In this case the current is determined by the 
probability for electrons to tunnel from the STM-tip into 
the CNT, and one could expect a characteristic length of 
about 0.1 nm. 

To conclude, we expect that vibrations of a suspended 
carbon nanotube can be self-excited in the DC-current bi- 
ased regime if a sufficiently large external magnetic field 
is applied perpendicular to the nanotube axis. We sup- 
port this claim with an analysis of a classical model of 
such a system and show that when the dissipation is suffi- 
ciently low, the static state of the nanotube becomes un- 
stable at reasonably weak (10-100 mT) magnetic fields. 
This magnetomotive instability develops into a steady 
state characterized by pronounced nanotube vibrations. 
We have also demonstrated that one effect of the nan- 
otube vibrations is to change the magnetic field depen- 
dence of the average voltage drop across the nanotube. 
This phenomena may be used for the experimental de- 
tection of the magnetomotive instability. 

Other types of active oscillators based on sus- 
tained self-oscillations of a on suspended CNT have 
been proposed theoretically^"^ and recently realized 
experimentally-^^. The approach in Ref. 24 relies on dis- 
tance dependent electron injection from a STM tip into 
a doubly-clamped CNT, while in Ref. '^"^ distance de- 
pendent field emission of electrons from a singly-clamped 
CNT to the electrode provides a feedback mechanism 
that leads to controllable mechanical self-oscillations 
from a single DC voltage supply. The drawback of both 
these devices is that they require very precise geometry 
control for their operation and that the short distance 
between the CNT and the controlHng electrode greatly 
limits the amplitude of vibration. The advantage of 
the device considered here is that the main element is 
a suspended-channel CNTFET — fabricated in numer- 
ous experimental laboratories — in which the CNT can 
be hundreds of nanometers from the gate electrode, and 
which requires lower voltages for operation. 
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